COSMIC RAYS IN THE ENERGY RANGE OF THE KNEE 
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A brief motivation for studying cosmic rays at energies 10^"* ^ E ^ 10^^ eV 
is given. Besides astrophysical interests in identifying and understanding their 
sources, there are also particle physics aspects related to their transport properties 
in the galaxy or to their detection via extensive air showers. The KASCADE 
air shower experiment taking data at Forschungszentrum Karlsruhe (Germany) 
provides important information to both of these topics of particle-astrophysics 
research. A target of particular interest is the so-called "knee" in the cosmic ray 
energy spectrum at E]^ ^ 4 ■ 10^^ eV. Recent results adding knowledge to an 
understanding its origin will be presented as well as results on tests of high-energy 
hadronic interaction models required for air shower simulations. 



1 Introduction 

The origin and acceleration mechanism of ultra-high energy cosmic rays have 
been subject to debate for several decades. Mainly for reasons of the power re- 
quired to maintain the observed cosmic ray energy density of Ecr ~ 1 eV/cm^ 
the dominant acceleration sites are generally believed to be supernova rem- 
nants (SNR). Charged particles mainly originating from the surrounding in- 
terstellar medium of the pre-supernova star may get trapped at the highly 
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supersonic shock wave generated by the SN explosion. Repeatedly reflec- 
tions on both sides of the shock front lead to an acceleration by the so-called 
'first order Fermi mechanism'. Naturally, this leads to a power law spectrum 
dj/dE (X E^'' as is observed experimentally. Simple dimensional estimates 
show that this process is limited to i?max ^ Z x {p x B), with Z being the 
atomic number of the cosmic ray (CR) isotope and p, B the size and mag- 
netic field strength of the acceleration region. A more detailed examination 
of the astrophysicaLparameters suggests an upper limit of acceleration of 
-Emax ^ Z X 10^^ eVUfl Curiously, the CR spectrum steepens from 7 ~ 2.75 
to ~ 3.1 at £' ~ 4 X 10^^ eV which is called the 'knee'. The coincidence thus 
may indicate that the 'knee' is related to the upper limit of acceleration. 

An alternative interpretation of the knee is that it may reflect a change in 
the propagation of CRs from their sources to the solar system. Propagation 
effects are studied mostly by balloon or space borne experiments measuring 
abundances of secondary to primary CR isotopes at energies up to some GeV. 
A prominent example is the B/C-ratio. Analysing such data yields the 'escape 
time' of CRs from our galaxy which scales with energy according to Tesc oc E~^ 
with 6 « 0.6. A simple extrapolation of Tesc to higher energies breaks down 
at iJ « 3 • ip^^ eV, because cTosc ^ 300 pc which is the thickness of the 
galactic diskD. The value corresponds to just one crossing of the disk and 
would give rise to significant anisotropics with respect to the galactic plane 
when approaching this value. Similarly as to the process of acceleration at 
SNR shocks, the process of galactic containment is closely related to magnetic 
field confinement, i.e. in addition to anisotropics one again expects E^^^ oc Z. 

A picture related-to both of these interpretations has been proposed by Er- 
lykin and Wolfendal^. They consider the knee as a superposition of a weakly 
energy dependent galactic modulation with additional prominent structures 
in the flux spectrum caused by a single near-by object. This so-called 'sin- 
gle source model' assumes that a shock wave of a recent nearby supernova 
which exploded some 10,000 years ago at a distance of a few hundred parsecs 
currently propagates, or has recently propagated through the solar system 
causing distinct peaks of elemental groups in the energy spectrum. However, 
there is some controversy whether or not the statistics of presently existing 
data gives support to the model. 

Recently, a more particle pkysics motivated picture of explaining the knee 
was put forward by WigmanaJ. He suggested that inverse beta decay of 
protons with massive relic neutrinos according to p + — > n + could 
destroy protons. Simple kinematics shows that this channel is open for Ep > 
1.7- 10-^^eV/my(eV). Thus, the knee energy E^kncc — 4 PeV would correspond 
to an electron neutrino mass of 0.4 eV, a value presently not excluded by any 
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other observation or experiment. However, 'eating' sufficiently large amounts 
of protons by such a process requires extraordinary high local densities of relic 
neutrinos, even if possible gravitational trapping is considered. 

It has also been suggested that the knee is not a property of the primary 
energy spectrum itself, but rather may be caused by changing high-energy 
interactions in the earth atmosphereH. Producing a new type of a heavy par- 
ticle in the first interactions escaping unseen by air shower experiments could, 
in principle, mimic a break in the spectrum. From the particle physics point 
of view this is not completely ruled out as the centre-of-mass energy avail- 
able at the knee is above Tevatron energies. Different from the astrophysical 
interpretations of the knee, one would now expect to see the knee energy of dif- 
ferent primary particles being displaced by their mass number A rather than 
by their charge Z. This is understood from the nuclear reaction mechanism 
being governed by the energy per nucleon E/A of the incident particle. 

A key observable for understanding the origin of the knee and distinguish- 
ing the SN acceleration model from the other proposed mechanisms, is thus 
given by CR energy spectra of different elemental groups or, if measured more 
inclusively, by the mass composition of CRs measured across the knee (see e.g. 
Ref.El). This is to be complemented by measurements of anisotropics in the 
arrival directions of CRs. Finally, experiments should also aim at testing 
hadronic interaction models by means of air shower data. 

At present, beyond the knee little is known about CRs other than their all- 
particle energy spectrumEI. The low flux of particles 1 m~^-year~^ above 
the knee) puts strong demands on the collection power of the experiments, 
such as can at present only be achieved by extended air shower (EAS) arrays 
at ground level. The development of an EAS is a process of truly multi- 
particle dynamics with frequent high-energy hadronic and electromagnetic 
interactions taking place. At primary energies of lO""^^ eV about 10^ particles 
reach sea-level. These are mostly photons and electrons plus some fraction of 
high-energy muons and hadrons, spread out over an area of several hectares. 
Sampling detector systems with typical coverages of less than one percent can 
thus be used for registration of such EAS. However, this indirect method of 
detection bears a number of serious difficulties in the interpretation of the 
data and requires-detailed modelling of the air shower development (see talk 



caused by the employed high-energy interaction models and by inevitable 
EAS fluctuations need to be considered. The basic concept of the KASCADE 
experiment is to measure a large number of EAS parameters in each individual 
event in order to verify the consistency of EAS simulations and to determine 
the energy and mass of the primary particles on a reliable basis. It is hoped 
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Particularly, systematic effects 




Figure 1. Schematic layout of the KASCADE experiment (left), with its streamer tube 
tracking system (top right) and central detector (bottom right)E|. 



that the new data will substantially improve our knowledge of the origin of 
high-energy cosmic rays. 

2 KASCADE Experiment 

KASCADE (Karlsruhe Shower Core and Array Detector) is located at the 
laboratory site of Forschungszentrum Karlsruhe, Germany (at 8° E, 49° N, 
110 m a.s.l.). In brief, it consists of three major components (see Fig.|l|); 

1. A scintillator array comprising 252 detector stations of electron and muon 
counters arranged on a grid of 200 x 200 m^ and providing in total about 
500 m^ of e/7- and 620 m^ of /i-detector coverage. The detection thresh- 
olds for vertical incidence are Ee > 5 MeV and iJ^t > 230 MeV. 

2. A central detector system (320 m^) consisting of a highly-segmented 
hadronic calorimeter read out by 40,000 channels of warm liquid ion- 
ization chambers, a layer of scintillation counters above the shielding, a 
trigger plane of scintillation counters in the third layer and, at the very 
bottom, 2 layers of positional sensitive MWPC's, and a streamer tube 
layer with pad read-out for muon tracking at i?^ > 2.4 GeV. 
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3. A 48 X 5.4 tunnel housing three horizontal and a two vertical layers of 
positional sensitive limited streamer tubes for muon tracking at E^^ > 0.8 
GeV. 

More details about the experiment can be found in Refs. illl. First cor- 
related data have been taken with some parts of the experiment in 1996 and 
with its full set-up since 2000. At present, more than 300 Mio. events have 
been collected in a very stable mode and with a trigger threshold of the array 
corresponding to -E ~ 4 • 10^"* eV. 

3 Tests of High-Energy Hadronic Interaction Models 

The observation of EAS provides an opportunity to study global proper- 
ties of hadronic interactions in an energy range not accessible to man-made 
accelerators. For example, the centre-of-mass energy at the Tevatron col- 
lider corresponds to a fixed target energy in the nucleon-nucleon system of 
Ep ~ 1.7 • 10^^ eV. Even more importantly, the diffractive particle produc- 
tion dominating the energy flux in the forward region and influencing the 
EAS development most strongly, was^udied experimentally only at compar- 
atively low energies of ~ 10 GeVB. Most of the beam energy in present 
collider experiments, however, remains unseen. For example, the UA5 ex- 
periment could register up to 30% of the total collision energy at y/s = 0.9 
TeV, while the CDF detector registers only about 5% at ^/s — 1.8 TeV. 
Hence, hadronic interaction models applied to higher energies and to parti- 
cle production in the forward region rely on extrapolations and may cause 
systematic uncertainties in simulations of EAS. Additional uncertainties arise 
from simulations of p-nucleus and nucleus- nucleus collisions including a pos- 
sible formation of a quark-gluon plasma. Again, such data are important for 
EAS interpretations but have been studied only at low energies in the past 
(SPS and ISR at CERN). Only very recently, RHIC data at = 200 GeV 
have become available and will be very helpful in this respect. Apart from 
different theoretical extrapolations, an additional source of uncertainty for all 
models originates_from the experimentally deduced energy dependent inelas- 
tic cross-sectionsEn At Tevatron energies, the proton- antiproton cross section 
is not known to better than 5 % at best. As will be shown below, these un- 
certainties are amplified when predicting absolute hadron fluxes at sea-level 
by means of EAS simulations. Thus, one may take advantage of this strong 
dependence and perform stringent tests of models on the basis of EAS data. 
The idea is to use data of absolute CR fluxes up to several TeV of energy, as 
obtained from balloon and satellite experiments at the top of the atmosphere, 
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Figure 2. Trigger rate vs hadron rate in KASCADE. The top panel compares different 
hadronic interaction models to the experimental data. The systematic uncertainty, mostly 
given by the absolute flux uncertainty of the direct experiments, i34jidicated by the dotted 
line. The lower panel shows results obtained fre«n|-^G QGSJETllj model with modified 
inelastic cross section and diffraction dissociationtjo. 



and propagate these particles (taking into account their energy distribution 
and chemical composition) through the atmosphere employing the CORSIKA 
simulation packagell3. At the level of the KASCADE experiment, we then 
ask for triggers released by either high energy hadrons {Eh > 90 GeV) or a 
minimum number of 9 muons detected in the central detector of KASCADE. 
These simulated inclusive trigger rates, employing different hadjixmic interac- 
tion models, are then compared to actual experimental datatJtj. Figure || 
shows the results. None of the predictions agrees well with the experimental 
data, particularly the hadron rates are overestimated by up to a factor of 3. 
Furthermore, there are also large differences found between the models. This 
convincingly proves the sensitivity of the experimental observable to details 
of the interaction models. 

Earlier investigations of high-energy hadrons observed in the shower core 



ISMD2000: submitted to World Scientific on February 1, 2008 6 



of the KASCADE calorimeteiO lead to the conclusion that QGSJET pro- 
vides the best overall prescription of EAS data. This was confirmed alsn 
by an independent analysis including data from different EAS experimentstS. 
Therefore, several modifications were applied to the QGSJET model in order 
to study their influence to the predicted trigger rates. Results are presented 
in the lower panel of Fig.||. Increasing the inelastic proton-air cross section 
by 5 and 10 % reduces the predicted hadron rate by approx. 27 % and 54 %, 
respectively. Similar effects are obtained by lowering the diffraction dissoci- 
ation by up to 10 % of the inelastic cross section. Clearly, there is a strong 
sensitivity to these parameters. Since the total inelastic cross section appears 
to be the better known quantity, we conclude that the diffraction dissociation 
is overestimated in the simulated hadron-nucleus interactions by about 5 % of 
the inelastic cross section (Tinei(p-air). Further studies with new models and 
with high energy hadrons observed in the core of EAS are in progress and will, 
hopefully, help to further improve our understanding of hadronic interaction 
processes in EAS. This will be important also for interpretations of EAS data 
in terms of primary energy and chemical composition. 

4 Energy Spectrum and Chemical Composition 

Extracting the primary energy and mass from EAS data is not straightfor- 
ward and a model must be adopted to relate the experimental observables to 
properties of the primary particle. As shall be discussed below, the analy- 
sis is complicated by large fluctuations of EAS observables and by the fact, 
that virtually all of these observables are sensitive to changes in the mass and 
energy of the primary particle. 

Air shower measurements generally include particle measurements at 
ground and are_sometimes complemented by detection of Cherenkov- and flu- 
orescence lighto. Primary observables of KASCADE are lateral particle den- 
sity distributions, Pe,pi,h{r), and their integrals Ne,^j.,h = Jl'^ 2'Krpe,^j.,h{r)dr, 
yielding total (if extrapolations ri — > and r2 ^ oo are made) or truncated 
particle shower sizes. Truncated shower-sizes, N^^^y^ , have been introduced 
by the KASCADE collaboration in order to avoid systematic uncertainties re- 
sulting from extrapolations of Pe.fi.hir) to large distances not covered by the 
experiment. Reconstruction techniques of the individual particle densities, 
shower core position, and shower direction have been discussed in Ref. E2l. As 
simulations predict and as it is obvious from very basic considerations, the 
weighted sum of particle shower sizes provides an estimate of the primary 
energy, i.e. E ^ a - logNjf + b ■ logN^'^, with a and b determined from sim- 
ulations. On the other hand, the ratio Yiatio — \og / log N^'^ provides a 
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reliable estimate of the primary 
masfcil. This is mostly because 
showers of heavy primaries will, 
on average, develop higher in the 
atmosphere so that the electro- 
magnetic component suffers more 
absorption in the atmosphere as 
compared to a shower induced by 
a primary protonE2l. An example 
is presented in Fig.||. Clearly, a 
superposition of different particle 
types (represented by the hues) 
is required to account for the ex- 
perimental distribution. Notewor- 
thy, the left and right hand tails 
of the experimental distribution 
are nicely described by the pro- 
ton and iron simulations, respec- 
tively. Performing such an analy- 
sis for different bins of shower sizes 
yields an increasingly heavier composition for energies above the knee". 

The separation of primary particle types by Yratio may also be used to 
disentangle the primary flux spectrum into 'light' and 'heavy' components. 
An example is shown in Fig.^. Here, the distribution of muon densities mea- 
sured in the trigger plane of the central detector at a distance of 45.5 m from 
the shower core is plotted for all showers and for those where Fratio < 0.75 
('electron-rich') and Fratio > 0.75 ('electron-poor'). Interestingly, the knee 
structure is seen for jmselected and electron-rich showers but is absent for 
electron-poor showersE3c3. Realizing that the muon density at given distance 
to the shower core is linked to the primary energy, the figure proves in a very 
model independent way that the knee structure appears to be a feature caused 
by the light component only. Again, we like to stress that is directly mea- 
sured by the experiment and that any meaningful model will predict electron- 
rich and -poor showers for light and heavy primaries, respectively. 

This finding is confirmed by combined energy and composition analyses 
of the shower size spectra of KASCADE, albeit then based on quantitative 
Monte Carlo simulations. A compilation of the all-particle energy spectrum 
reconstructed from various experiments is presented in figure I The agree- 
ment appears reasonable and deviations are mostly explained by uncertainties 
in the energy scale by up to 25%, e.g. CASA MIA datao were shifted up- 
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Figure 3. Example of an event-by-event distri- 
bution of the muon/electron ratio measured by 
KASCADE for a primary energy of E ~ 310^^ 
eV. The lines show results of CORSIKA simu- 
lations employing the QGSJET model and are 
performed for different primary particles. 
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Figure 4. Distribution of muon densities (_Ep > 490 MeV) measured in the trigger plane 
at a distance of 45.5 m from the shower axis. The 'all'-particlc spectrum is shown as well 
as spectra forr^lectron-rich' and '-poor' showers. The black triangle indicates the position 
of the kneeB'El. 
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Figure 5. Primary energy spectrum as obtained from different experiments. The lines 
represent KASCADE results based on a simultaneous fit to electron and muon shower sipep 
assuming the all-particle spectrum to be a superposition of a proton and iron component^. 



wards in energy by 20% to yield a better agreement to the other data-sets. 
The shift is hkely to be explained by the interaction model SIBYLL 1.6cj em- 
ployed by the authors of Ref. c3 but which_has been proven to provide only a 
poor description of the experimental datallZl. The lines correspond to a simul- 
taneous fit of the electron and muon size spectra of KASCADE, assuming the 
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all-particle spectrum to be described by a sum of proton and iron primarie: 




Energy and mass dependent fluctuations of the electron and muon shower size 
are taken into account in this analysis and, in fact, turn out to be important 
for any quantitative discussion. Again, the knee is only reconstructed for the 
light component and no indication of a break is seen in the heavy component 
up to about 10^'' eV. Obviously, the different spectral shapes of the 'light' and 
'heavy' induced air showers will result in an increasingly heavier composition 
above the knee, such as is found also by other experiments. The analysis 
is an ongoing process with improved statistics and reconstruction techniques 
making use of the multi-detector capabilities of KASCADE.EZI 



5 Conclusions and Outlook 

KASCADE data provide mounting evidence for the knee being caused by the 
light particle component only. Up to an energy of ~ 10^^ eV, which is the 
upper limit of the present KASCADE experiment, no indication of a similar 
structure is found in the iron-like distribution. However, comparison to the 
all-particle spectrum deduced from other experiments, suggests the necessity 
of a similar break in the heavy component just above 10^^ eV. This important 
finding gives direct support to the picture of acceleration or galactic confine- 
ment in magnetic fields. To discriminate the two pictures, high statistics data 
on anisotropy will be required and are currently being analysed. Furthermore, 
an extension of the mass-selective measurements to higher energies appears 
mandatory in order to test the hypothesis of a 'Fe-knee' in more detail. This 
should go in parallel with further tests of hadronic interaction models in order 
to improve the reliability of the EAS-data analysis. 

The KASCADE and EAS-TOPEa collaborations have just started a com- 
mon effort to realise these goals already in the near future. The EAS-TOP 
scintillators are currently positioned next to the KASCADE experiment at 
the site of Forschungszentrum Karlsruhe providing a 12 times larger accep- 
tance as compared to the present set-up. Generating a common trigger for 
the combined experimental installations will allow taking advantage of the full 
multi-detector capabilities including hadronic measurements as well as muon 
counting, tracking, and timing capabilities. The new installation is foreseen 
to start taking data in spring 2002. It is expected that the upcoming KAS- 
CADE and future KASCADE-GRANDE data will allow extending the tests 
of hadronic interaction models to higher energies and to substantially improve 
our understanding of the origin of the knee in particular and about the origin 
and acceleration of cosmic rays in general. 
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